Central depression in nuclear density and its consequences for the shell structure of 

superheavy nuclei. 
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The influence of the central depression in the density distribution of spherical superheavy nuclei 
on the shell structure is studied within the relativistic mean field theory. Large depression leads 
to the shell gaps at the proton Z = 120 and neutron iV = 172 numbers, while flatter density 
distribution favors N = 184 for neutrons and leads to the appearance of a Z = 126 shell gap and 
to the decrease of the size of the Z = 120 shell gap. The correlations between the magic shell gaps 
and the magnitude of central depression are discussed for relativistic and non-relativistic mean field 
theories. 
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The question of the possible existence of shell- 
stabilized superheavy nuclei and the precise location of 
the magic spherical superheavy nuclei has been in the 
focus of the nuclear physics community for more than 
three decades Unfortunately, the various theoretical 
models do not agree with respect to the magic shell gaps 
in superheavy nuclei. The proton numbers Z = 114, 120 
and 126 and the neutron numbers N = 172 and 184 are 
predicted by different models and parametrizations 0. 
The same models reproduce the known magic numbers 
in lighter systems. The predicted magic numbers are of 
decisive importance for the experimental search of su- 
perheavy nuclei. In such a situation it is necessary to 
understand what makes the predictions of the models so 
different in the region of superheavy nuclei. One of the 
reasons is the appearance of a central depression in the 
nuclear density II, 3, which is studied in this paper. 

The first predictions of superheavy nuclei were based 
on the shell correction method, which assumes a single- 
particle potential with a flat bottom. Nowadays, these 
calculations predict Z = 114 and N=184. The micro- 
scopic self-consistent models start either from effective 
nucleon-nucleon interaction (models based on the Skyrme 
and Gogny forces) or effective exchange of mesons by nu- 
cleons (relativistic mean field (RMF) theory). Although 
based on more fundamental principles, these models do 
not agree among each other in predicting the magic shell 
gaps of superheavy nuclei. In part, this is related to 
the fact that the reliability of different parameterizations 
of these models is verified only by comparing theoreti- 
cal and experimental binding energies and its derivatives 
(separation energies, the S2 n ,2p(Z, N) quantities) and de- 
formation properties These observables are not very 
sensitive to the energies of the single-particle states. For 
example, it was shown in Ref. [5| that the NLSH and 
NL-RA1 parameterizations of the RMF theory provide 
a reasonable description of these quantities in deformed 
actinide region despite the fact that the single-particle 
energies are poorly reproduced. Accurate single-particle 



energies are crucial for predicting the shell gaps in su- 
perheavy nuclei. However, the accuracy of the descrip- 
tion of the single-particle states in deformed region of 
the heaviest actinides has been tested only for few pa- 
rameterizations of the RMF theory || and Skyrme SLy4 
functional [f|. The shell correction approach is based on 
phenomenological potentials which reproduce the single 
particle levels of the actinides best. However, the as- 
sumption of a flat bottom radial profile is a severe source 
of error when extrapolating to spherical superheavy nu- 
clei. 

Self-consistent microscopic calculations find a central 
depression in the nuclear density distribution 0,0], which 
generates a wine-bottle nucleonic potential. Its magic 
numbers differ from the ones of the phenomenological 
flat-bottom potentials. The present manuscript studies 
the influence of this depression on the shell structure of 
spherical superheavy nuclei. As a theoretical tool we use 
the RMF theory for spherical nuclei without pairing |7| 
and the relativistic Hartree-Bogoliubov (RHB) theory [8[. 

Fig. Q] compares the single-particle spectrum of a wine- 
bottle potential (g-s) with the one of a flat-bottom poten- 
tial (exc-s). The details of the potentials will be discussed 
later. The differences are easy to understand. The high- 
j orbitals are localized mostly near the surface, whereas 
the low-j orbitals have a more central localization. As 
compared to a flat-bottom potential, in general the high- 
j orbitals are more and the low-j orbitals are less bound 
in an attractive wine-bottle potential. In the following 
it is useful to distinguish between the groups of low- 
j and high-j single-particle states. Filling up a low-j 
group with nucleons increases the density near the cen- 
ter, whereas filling a high-j group increases the density 
near the surface. As we shall demonstrate, the occupa- 
tion of these groups determines the radial profile of the 
neutron and proton densities and potentials. 

We start with 208 Pb. The RMF theory provides a good 
description of the experimental charge density distribu- 
tion of this nucleus [{J, • With increasing neutron and 
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FIG. 1: Single-particle spectra of ground state (indicated 
as 'g-s') and excited (indicated as 'exc-s') configurations in 
292 120i72 system obtained in the RMF calculations with the 
NL3 force. Solid and dashed lines are used for positive and 
negative parity, respectively. Solid and open circles indicate 
the occupied and empty orbitals, respectively. In the ground 
state, all subshells below Z — 120 and N = 184 are fully 
occupied. For the g-s configuration, the spin-orbit partners 
3pi/2 5 3^3/2 and 3d 3 / 2 , 3d 5 / 2 show up at very close energy. In 
the excited configuration, only 12 particles are excited from 
the subshell ^lj'15/2: 4 particles still reside in this subshell. 
The spherical shell gaps of interest are indicated. 



proton numbers the corresponding densities are modified 
in the way shown in Fig. [51 Between Z — 82 and Z = 106 
the protons fill the high-j group 7rlz 13 / 2 , 7r2/ig/2, between 
Z = 106 and Z = 120 they fill the medium-j group 
^2/7 /j) 7i"2/ 5 / 2 , and between Z = 120 and Z = 126 they 
fill the low-j group 7r3p 1 / 2 , n3p3/2- The variations of the 
proton density are seen in the N = 172 isotones. The 
filling of the high-j group 7rli 13 / 2 , 7r2/i g / 2 increases the 
density at the surface (compare Z = 82 and Z = 106 
in Fig. [5]). The filling of medium-j group 7r2/ 7 / 2 , 7t2/ 5 / 2 
adds to the density between central and surface areas 
(see Z = 120). Finally, the filling of the low-j group 
^Pi/2j 7r ^P3/2 adds to the density in the near-central re- 
gion of nucleus (see Z = 126). 

The analogous grouping into low-j and high-j sub- 
shells in the neutron subsystem is illustrated in Fig. [3] 
The variation of the neutron density generated by filling 
these groups is seen most clearly in the Z = 106 iso- 
topes. Filling the high-j group vliu/2, fl.7'15/2, v2g$/2 
increases the density near the surface. Filling the low- 
j group v3d 5 / 2 , vZd%/2, visi/2 increases the central den- 
sity and filling the high-j group vlj 15/2 , ^2/i 11/2 , vlkn/ 2 
adds to the surface region again. Analyzing the 
published results, we found that the grouping into 
high/medium/low- j subshells shown in Fig. [5] appears 
in all models/parametrizations (cf. Fig. [1] in the present 



manuscript and Figs. 4, 9, 13, and 15 in Ref. @). 

As seen in Fig.[3J the combined occupation of the high- 
j neutron subshells 2g g / 2 , lji5/ 2 , H11/2 [and medium-j 
2<77/ 2 ] and proton l/ig/ 2 and Hi 3 / 2 [and medium-j 2/7/2] 
subshells leads to a central depression in the nuclear 
density between Z = 106 and Z = 120 and N = 164 
and N = 172, which is especially pronounced in the 
Z = 120, N = 172 system. As seen from the density 
variations in Fig. [51 the proton subsystem plays a larger 
role in the creation of the central depression. This result 
differs from the results of the Skyrme calculations with 
the SkI3 parametrization [3j , the authors of which claim 
that the central depression is mainly due to the occupa- 
tion of neutron subshells. The appearance of the central 
depression is a consequence of the different density dis- 
tributions of the single-particle states: high-j orbits are 
located near the surface and low-j orbits near the center. 
This generic feature is dictated by the nodal structure of 
the wave functions in a leptodermic potential. Hence, the 
high-j proton and neutron orbitals will modify the radial 
profile in a comparable way. However, the high-j proton 
orbitals should be more efficient, because the Coulomb 
interaction pushes them to larger radii. 

Let us study the interplay between the geometry of 
the single-particle orbitals, the appearance of the central 
depression in the density and the shell structure in more 
detail. One possibility is to generate a flatter density dis- 
tribution in the central part of nucleus by exciting parti- 
cles from high-j subshells to low-j subshells. Fig.[T]shows 
an example of such an excitation in the 292 120i7 2 system. 
Here 12 neutrons are excited from the ljis/ 2 subshell into 
3d 5 /2, 3e? 3 / 2 and 4s 1 / 2 subshells. In this excited (called 
'exc-s') configuration the neutron density distribution in 
the central part of nucleus is much flatter than in the 
ground state (called 'g-s') configuration, and its profile is 
very similar to the one in 208 Pb (Fig.^j). The changes in 
neutron density are fed back to proton density, because 
the isovector interaction tries to keep them alike. As 
a consequence, the proton density distribution becomes 
also flatter (Fig. [3^,) but density fluctuations due to shell 
effects remain visible. The nucleonic V + S (Fig. [3fc,d) 
potential and the (V — S) (Fig. [3b, f) potential, which in 
first approximation is related to the spin-orbit potential 
via V ls (r) = ^(V(r) - S(r)) 0, reflect the density 
change: they becomes flatter in the central part of nu- 
cleus. This effect is especially pronounced in the proton 
nucleonic potential: the 'wine-bottle' radial shape is re- 
placed by a 'flat-bottom' one. Another consequence of 
this excitation is an increase of the surface diffuseness 
both in the densities and in the potentials. 

The various RMF forces are characterized by different 
compression moduli (NL-Z [K^ = 173 MeV], NL3 
[Koo = 272 MeV] and NLSH [K^ = 355 MeV]). As 
expected, the magnitude of central depression in densities 
and potentials increases with the decrease of compression 
modulus. However, the changes in the energies of single- 
particle states, densities and potentials induced by our 
probing particle-hole excitation do not depend sensitively 
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FIG. 2: The evolution of proton and neutron densities with the changes of proton and neutron numbers. Arrows indicate the 
group of single-particle subshells which become occupied with the change of the nucleon number. The figure is based on the 
results of spherical RMF calculations without pairing employing the NL3 parametrization. The shaded background is used for 
nuclei located beyond the proton-drip line. If the indicated configuration is not lowest in energy, its excitation energy (in MeV) 
is given by E*. 



on the compressibility. 

As a result of the flattening of the nucleonic potential 
the energies of the single-particle states are changed as 
described above (see Fig. [IJ . The shifts are larger in the 
proton subsystem because the proton central potential is 
more flattened than the neutron one. We see the Z = 126 
proton gap emerging and the size of the Z = 120 gap 
decreasing. To a lesser extend, the N = 172 neutron gap 
decreases and the N = 184 gap increases. The flattening 
of the (V — S) potential increases the splitting of the spin- 
orbit pairs [7r3p 1/2 , 7r3p 3/2 ], [f3d 3 / 2 , vM h / 2 ], and [7r2/ 5/2 
and 772/7/2]. The spin-orbit splitting of the last pair of 
orbitals generates the Z = 114 shell gap predicted by a 
number of models. The present results clearly show that 
a flatter density distribution leads to a larger splitting 
between these orbitals. 

We have studied further excitations that induce a flat- 
ter density distribution. In all cases we found the above- 
mentioned dependence of the size of the Z — 120, 126 and 
N = 172, 184 shell gaps on the magnitude of the central 
density depression. Our results are consistent with the 
HFB studies with the Gogny D1S force, which employ 
an external potential in order to induce the central de- 
pression [ij: large N = 184 and Z — 126 shell gaps were 
found for the values of the external potential that gen- 



erate a fiat density distribution and large Z — 120 and 
N = 172 shell gaps for the values that generate a central 
depression (see Fig. 2 in Ref. [J]). 

Due to the isovector force, which tries to keep the neu- 
tron and proton density profiles alike, there is a mu- 
tual enhancement of the Z = 120 and A" = 172 gaps, 
both being favored by the wine-bottle potential, and of 
the Z = 126 and A" = 184 gaps, both favoured by the 
flat bottom potential. For the same reason the gaps are 
smaller for the combination Z — 126 and A~ = 172, and 
the Z = 120 gap does not develop for A" = 184. This 
behavior is not expected to depend much on the den- 
sity functional chosen. Indeed, a number of Skyrme cal- 
culations (SkI3, SkI4, Skll, SLy6), which show a large 
Z = 120 gap in the 292 120u 2 system, do not show dou- 
ble shell closure at Z = 120, N = 184 13]. These generic 
features are also seen in the calculations with Gogny D1S 
force (Fig. 2 of Ref. [HI), with the Skll (Fig. 2 in Ref. 
0), and SkI3, SkI4, SkP (Figs. 6, 7, and 8 in Ref. g) 
Skyrme forces, and with the RMF NL3 and NL-Z2 forces 
(Fig. 2 in Ref. 0). 

Let us consider the 292 120u 2 system, which is a dou- 
bly magic superheavy nucleus in RMF theory. Both 
relativistic and non-relativistic (Gogny D1S, Skyrme 
parametrizations with low isoscalar effective mass m*/m 



4 



0.1 



CD 



proton 




neutron 


1 1 1 1 

1 (a) J 




1 1 1 

""N? (b) - 

\\ 


/ ■ 

/ / v.- 




\\ J 


: /V- 208 Pb : 




\\ 

\\ ~ 
i . -\l 




Radius r [fm] 

FIG. 3: Density distributions (upper row), nucleonic (V + S, 
middle row) and V — S (bottom row) potentials in the ground 
(g-s) and excited (exc-s) configurations of 292 120i72 . S and V 
are attractive scalar and repulsive vector potentials, respec- 
tively. The left column shows the proton system and the right 
the neutron system. The proton and neutron density distri- 
butions of 208 Pb are shown in upper panels for comparison. 
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FIG. 4: Neutron and proton densities of 292 120i72 obtained 
in different models/parametrizations. The densities obtained 
with Skyrme and Gogny forces are taken from Refs. pH Hl|. 



such as SkI3, SLy6 [3]) show a pronounced central de- 
pression (see Fig. 2]). These density functional are 
characterized by similar values of m*/m (Gogny D1S 
[m*/m = 0.67], Skyrme SkI3 [m*/m = 0.57] and SLy6 
in */m = 0.69] 0. These values should be compared 
with RMF Lorentz effective mass of the nucleon at the 
Fermi surface m*{kp)/m w 0.66 [|[, since effective mass 
is momentum-dependent in the RMF theory (la |. The 
central depression is much smaller in the Skyrme calcu- 
lations with SkP (Fig. SI) and SkM* forces (Ref. 0) which 
have high values of isoscalar effective mass m* jra — 1 and 
m* jm = 0.789, respectively. The development of a more 
pronounced central depression for the density function- 
als with low effective mass may be understood as follows. 
In the surface region, the ratio m*/m changes from its 
value < I in the interior to 1 in the exterior. Classi- 
cally, nucleons with given kinetic energy are more likely 
to be found in regions with high effective mass than in 
the regions with low one because they travel with lower 
speed. This is reflected by the Thomas-Fermi expression 
for the nucleonic density p oc [2m*(eF — V)] 3 ^ 2 . The in- 
crease of the effective mass in the surface region favours 



the transfer of mass from the center there, which makes 
the above discussed polarization mechanism of the high- j 
orbitals more effective for functionals with low effective 
mass. Based on this argument we suggest that a flatter 
radial profile is a generic feature of the density functionals 
with an effective mass close to one. It would be interest- 
ing to investigate if the Skyrme functionals of this type 
systematically give flatter density distributions than the 
ones with a low effective mass. 

All experimentally known nuclei with Z > 100 are ex- 
pected to be deformed [l], Gil- The deformation leads 
to a more equal distribution of the single-particle states 
emerging from the high-j and low-j spherical subshells 
(see, for example, the Nilsson diagrams in Figs. 3-4 in 
Ref. [13]) than for spherical shape. Thus, the density 
profile of deformed nuclei is relatively flat |l8[, strongly 
resembling the one used in phenomenological potentials. 
This together with a careful fit of the single-particle en- 
ergies to the deformed nuclei in heavy actinide region ex- 
plains the success of the shell correction method [19l [20j | . 
However, this method neglects the self-consistent rear- 
rangement of single-particle levels due to the appearance 
of a central depression in spherical superheavy nuclei. 
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Thus predictions of magic numbers for superheavy nuclei 
within the shell correction method should be considered 
with caution. 

We have deliberately excluded from our study the 
forces NLSH and NL-RAl [RMF] and SkI4 [Skyrme], 
which give a Z = 114 shell gap in self-consistent cal- 
culations. This is because they provide a poor descrip- 
tion of either the energies of single-particle states in de- 
formed actinide A ~ 250 nuclei [5|, [2l| or of the spin- 
orbit splitting [3]. The energy splitting between de- 
formed states emerging from the 2/7/2 and 2/5/2 sub- 
shells is well described by the RMF NL1 and NL3 @ 
and Skyrme SLy4 Q forces, which give a small Z = 114 
shell gap. In our opinion, these results make the pre- 
dicted shell gaps at Z = 120, 126 and N = 172, 184 most 
likely. As discussed above, their appearance and combi- 
nation depends on the magnitude of the central depres- 
sion. The RMF theory gives a pronounced double shell 
closure at Z = 120, N = 172. The non-relativistic theo- 
ries (Gogny [ll[, Skyrme [IHQjj]) gi ye a l ar g e shell gap 
at N = 184 and less pronounced gaps at Z = 120 and 
126, the size of which strongly depends on neutron num- 
ber. For example, the Skyrme forces with high effective 
mass (SkM*, SkP) tend to predict a double shell closure 
at Z = 126, N — 184, while those with low effective mass 
(Skll, SkI3, SkM, SLy6) show a large gap at Z = 120 
for N — 172, which becomes smaller or disappears when 
approaching N = 184. 

In RMF theory, the N = 172 gap lies between the sub- 
shells v3d 5 / 2 and ^2<7 7 / 2 , which form a pseudospin dou- 
blet [HI]. The analysis of their deformed counterparts 
in the A ~ 250 region shows that the experimental en- 
ergy distance between the pseudospin partners i/l/2[620] 
and ^3/2 [622] is well reproduced, which supports the pre- 
dicted existence of a gap at N = 172 (see Fig. 28 in 
Ref. Q). However, taking into account the typical un- 



certainty of the description of the single-particle states in 
best-tested RMF parametrizations, one cannot exclude a 
large gap at N — 184 [f|. For this to take place, the 
energy of v^S\/2 state has to be overestimated by ap- 
proximately 1 MeV. 

In summary, the influence of the filling of the spherical 
subshells on the radial density profile and shell structure 
of superheavy nuclei has been studied. The occupation of 
high-j subshells decreases the density in the central part 
of the nucleus, the occupation of low-j subshells increases 
it. The polarization due to high-j orbitals generates a 
central depression of the density for nuclei with Z ss 120 
and/or N « 172, which is particularly pronounced for the 
combination Z = 120, N = 172, because both the proton 
and the neutron subsystems induce a central depression. 
This large central depression produces large shell gaps at 
Z = 120 and N = 172. The occupation of low-j orbitals 
by means of either multi-particle-hole excitations or of 
the increase of Z, N beyond Z = 120, N — 172 removes 
the central depression and reduces these shell gaps. The 
shell gaps at Z = 126 and N = 184 are favored by a 
flat density distribution in the central part of nucleus. 
The magnitude of central density depression correlates 
also with the effective mass of nucleons: low effective 
mass favors the large central depression. The similarties 
and differences between non-relativistic and relativistic 
mean field models in the predictions of magic shell gaps 
in spherical superheavy nuclei were discussed. 
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